Absolute Dynamical Limit to Cooling Weakly-Coupled Quantum Systems 
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Cooling of a quantum system is limited by the size of the control forces that are available (the 
"speed" of control). We consider the most general cooling process, albeit restricted to the regime 
in which the thermodynamics of the system is preserved (weak coupling). Within this regime, we 
further focus on the most useful control regime, in which a large cooling factor, and good ground- 
state cooling can be achieved. We present a control protocol for cooling, and give clear structural 
arguments, as well as strong numerical evidence, that this protocol is globally optimal. From this 
we obtain simple expressions for the limit to cooling that is imposed by the speed of control. 
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Preparing quantum systems in pure states is important 
for potential quantum technologies [TH1] , and this task is 
strongly linked to cooling to a (non-degenrate) ground- 
state: both require that all the entropy is extracted from 
the system. For this reason, there is presently a great 
deal of interest, experimental and theoretical, in cool- 
ing nano-mechanical resonators, and a number of cool- 
ing schemes of increasing effectiveness have been pro- 
posed [SHH] . Since the forces used to implement cool- 
ing are always limited, the question of the maximum 
achievable ground-state population for a given maximum 
control force is important both fundamentally and prac- 
tically. There are two structurally different regimes of 
cooling, that in which the dynamics of the thermal re- 
laxation is preserved under the control (weak coupling), 
and that in which it is changed (strong coupling). Here 
we address the problem of optimal cooling in the former. 

The complexity of the cooling problem — and the 
closely related state-preparation problem — is due to the 
interplay of both unitary and irreversible dynamics. It 
is not usually possible to prove the optimality of control 
protocols for complex dynamical systems, and this has 
hampered progress in determining the fundamental lim- 
its of quantum control. We adopt a heuristic approach, 
and explore whether we can use the structure of the cool- 
ing problem to guess an optimal cooling protocol under 
dynamical constraints. We test the optimality of this 
protocol by comparison with those found using numeri- 
cal optimization. We present very strong analytical and 
numerical evidence that we have explicitly constructed a 
globally optimal protocol, and thus determined the ab- 
solute dynamical limit to cooling in the dual regimes of 
weak coupling and good ground-state cooling. 

We consider the most general setting in which an TV- 
dimensional "target" system can be cooled: the target 
system is coupled to a second, M-dimcnsional "auxil- 
iary" system via an interaction Hamiltonian, Hi, whose 
eigenvalues we denote by HXj . This Hamiltonian, cou- 
pled with any trace-preserving operation on the auxiliary, 
implements the cooling process. The constraint we im- 



pose on the speed of control is that |Aj| < g, Vj, for some 
rate constant g. For a given experimental scenario, one 
calculates the eigenvalues Xj by i) determining the full 
Hamiltonian for the combined target and auxiliary sys- 
tems; ii) removing all matrix elements that act on only 
one of the subsystems; iii) calculating the eigenvalues of 
the matrix that remains. 

Our results apply to all cooling schemes, including 
those involving continuous measurement via a probe sys- 
tem [9]. However there is one class of schemes to which 
they are not immediately applicable: those in which the 
target is coupled directly to a near-continuum of states 
(e.g. a continuous measurement made directly on the 
target). This is because it is not yet clear whether our 
dynamical constraint is the right one from a practical 
point of view for this situation; the question of how it 
relates to the natural parameters for that case — the 
"damping rate" or "measurement strength" — is an in- 
teresting question for future work. 

To proceed we must chose some specific model of ther- 
malization (environmental noise), and the Redfield mas- 
ter equation is the obvious choice [10]: it gives an accu- 
rate description of any weakly damped Markovian quan- 
tum system, and weakly damped systems are the most 
important for quantum technologies. (The Redfield mas- 
ter equation reduces to the standard quantum optical 
master equations for the Harmonic oscillator and two- 
level systems [IS]-) We also start our analysis by allow- 
ing the auxiliary system to be ideal: it has large enough 
energy gaps that it has zero entropy at the ambient tem- 
perature (it is in its unique ground state), and is not itself 
subject to any damping or decoherence processes (we will 
see that if the auxiliary has zero entropy, damping it does 
not improve cooling). To a good approximation, present 
cooling methods for ions and nano-resonators are ideal 
in the first manner, but not the second. 

To cool the target, the control applied via the aux- 
iliary must take the target out of equilibrium with its 
bath, and the bath therefore induces an irreversible re- 
laxation of the system during the cooling process. We 
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FIG. 1. (a) Depiction of the density matrix for the joint state 
of the target and auxiliary. The auxiliary states are the "fast" 
index, illustrated here with M = 4. (b) A schematic depiction 
of the geometry of quantum dynamics, showing why rotations 
in the local subblocks do not change the angle-to-go for the 
cooling rotation. The "great circle" that cools is the dashed 
arrow, and the dashed ellipse is the motion of a local rotation. 

must understand how the control and damping act to- 
gether. Fortunately the regime in which cooling is most 
useful, and where present experiments operate, is that 
in which the target is cooled by a large factor / 3> 1, 
where / is the reduction in the average excitation num- 
ber. This requires that the damping (thermalization) 
rate of the system, 7, satisfies 7/ 'g(n + 1) <C 1 (where n 
is defined below). The means that we can analyze this 
regime by considering the dynamics to first order in our 
small parameter, and this will allow us to obtain ana- 
lytic results. If we define T>(c)p = {Sep + pSc)/2 — cpS 
for some operator c, then the master equation that de- 
scribes thermalization of a two-level system at tempera- 
ture T is given by p = — + n)T>(a) +nT>(a^)]p, where 
n = exp~ huJ / kT /(l — exp~ hw / kT ). Here hu is the energy 
gap, a is the lowering operator, and the equilibrium pop- 
ulation of the excited state is Pr = n/(l + 2n). For a 
harmonic oscillator the master equation is the same, but 
with the replacement a — > a, uj is now the frequency of 
the oscillator, and n now gives the average number of 
photons/phonons at equilibrium. 

Let us denote the energy levels of the target by |m), 
and those of the auxiliary by |j) x . If the auxiliary is 
prepared in the state |0) x , then the initial density ma- 
trix for the combined target and auxiliary is given by 
the matrix in Fig. [TJ Each subblock of this matrix is 
the full state-space of the auxiliary, and corresponds to 
a single state of the target. We see that the initial popu- 
lations each appear in the elements labelled by Ax, and 
are thus in different subblocks, and the task of cooling 
is that of transferring all the population to the subblock 
in the upper left-hand corner. The key observation we 
need is that cooling is a process of population transfer 
between orthogonal subspaces. Using the concept of ma- 
jorization, it is readily shown that the largest possible 
population in any subspace occurs in the basis in which 



the density matrix is diagonal [TT] , so we can expect that 
an optimal control protocol will keep the density matrix 
close to diagonal. 

Now that we have reduced the cooling problem to the 
transfer of population between orthogonal subspaces, in- 
sights from the geometry of quantum dynamics are now 
very useful [T2UT4] . The first of these is that given the 
constraint above, the fastest way to take any initial pure 
state to any other state is via a geodesic, the equivalent 
of a great circle in real vector spaces. The rotation angle 
along this great circle is determined by the inner product 
between the current and final states. The minimum time 
to get from any state |1) to an orthogonal state \2) is 
t = n/(2g), and is achieved by the simple Hamiltonian 
H = g(\l){2\ - |2)(1|). This tell us immediately that in 
the absence of damping (thermalization), the fastest way 
to perform the cooling operation is to rotate the states 
labelled by Aa; in Fig. [TJ respectively to each of the un- 
populated states labelled by Bx, at this rate. 

A second fact, illustrated in Fig. [ljb), is that local 
Hamiltonians of either subsystem cannot perform the 
population transfer needed to cool. Local Hamiltonians 
can only transform states to other states within specific 
subspaces - we will call these the local spaces of each sys- 
tem, and such population redistribution does not change 
the entropy of either system. Cooling therefore requires 
that we transfer population between two or more of these 
local subspaces (e.g between the A and B subspaces). 
Further, from the dynamical geometry, if Hi has par- 
tially rotated a state from one local space to another, no 
local Hamiltonian can change the angle (inner product) 
remaining between the rotated state and the destination 
B-state. Thus it cannot decrease the cooling time. 

Finally, by switching to the interaction picture, geome- 
try reveals immediately that any action of a local Hamil- 
tonian at any time can be undone, as far as the total pop- 
ulation transfer is concerned, by varying Hi with time: 
all we have to do is change the great circles upon which 
Hi rotates states, and this preserves the evolution of the 
total entropy of each system. Since we wish to obtain the 
optimal cooling over all possible time-dependent control 
protocols, this means that we can set the local Hamilto- 
nians to zero, with no loss of generality. With the above 
facts we can make the following general statement: 

Theorem: Trace-preserving operations on an auxiliary 
system can reduce, but not increase, the maximum pos- 
sible population transferrable to the ground state of the 
target system at any future time. 

Proof The population transferred to the ground state 
is the total population transferred between two orthog- 
onal local subspaces, and this transfer is achieved most 
rapidly by a set of 90° rotations on great circles connect- 
ing the two spaces. The action of a Hamiltonian local to 
the auxiliary cannot change the transfer rate. Further, i) 
any trace-preserving operation on the auxiliary is given 
by some unitary, U, acting on the auxiliary and an ar- 
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FIG. 2. (Color online) The performance of our conjectured 
optimal cooling protocol, compared against numerical opti- 
mization. The ground state population is Pg, and 1 — P g is 
shown vs. time for systems of size N, cooled by an auxiliary 
of size M, with control rate g — 1, damping rate 7 = 0.01, 
and initial thermal factors n. (a) Cooling using an ideal (un- 
damped) auxiliary, in which our protocol is optimal. Dark 
line: (N,M) = (2,3), n = 0.5; dashed line: (N,M) = (4,4), 
h — 0.5 light line: (N,M) = (4,4), n = 0.1; circles, squares, 
and diamonds: the corresponding results for numerical op- 
timization. Inset: numerical optimization beyond the opti- 
mal cooling time (see text), (b) Cooling the same systems 
as in (a), but the auxiliary is strongly damped with rate 
k = 1. Our protocol remains optimal for cooling a single 
qubit ((N,M) = (2,3)). 



bitrary third system, S3, and ii) as far as the target is 
concerned, U is local in the sense above. Thus U cannot 
increase the transfer rate, but it can reduce it: the inter- 
action between the target and auxiliary, Hi, can rotate 
on all great circles that connect local subspaces of their 
joint space, but it does not have the same access to all 
subspaces of the auxiliary and S3. After the action of U, 
Hi may no-longer be able to use the shortest paths for 
the transfer, and the transfer time will increase [T5]. □ 

Having established the Hamiltonian structure of the 
cooling problem, we now consider the continual thermal- 
ization of the target during the control process. We first 
present an argument that gives strong support for the fol- 
lowing simple and rather remarkable statement: all opti- 
mal cooling protocols will achieve the maximal ground- 
state population just prior to a time r = ir/ (2g), and will 
do so only when the system starts in equilibrium. This 
argument is as follows. The thermal master equation 
gives transition rates between the diagonal elements of p, 
and decay rates for all off-diagonal elements (in the en- 
ergy basis) . As soon as the ground-state population rises 
above the equilibrium value, there is a net thermal tran- 
sition rate out of the ground state to states orthogonal 
to it. The population taken out of the ground state can 



be returned by transferring it to the auxiliary, which, as 
established above, will require a minimum time r. Dur- 
ing this time, a further amount of population, P(t), will 
be taken out by thermalization. It is this amount of 
population that we will never be able to restrict to the 
ground state, since this population will again flow out 
by the time we have transferred it all back. Further- 
more, since population transfer is a rotation on the unit 
sphere, the rate of the increase of the ground-state pop- 
ulation goes to zero as t approaches r. Since the exit 
rate from the ground state is nonzero at t = r, no matter 
how large 3/7, there will always be a time slightly prior 
to r at which we lose by waiting longer. Finally, since 
the outward transition rate increases as the ground-state 
moves away from equilibrium, P(r) will only be minimal 
when the ground-state starts with its equilibrium value. 
This implies that the maximal cooling can only be ob- 
tained instantaneously; no steady-state cooling protocol 
can achieve this maximum. 

To arrive at our proposed optimal cooling Hamiltonian, 
we need to examine where thermalization places the pop- 
ulation that leaves the ground state. For two-level sys- 
tems and harmonic oscillators, thermalization involves 
transition rates only between adjacent energy levels, and 
so we focus on this here (the analysis generalizes easily). 
To first-order in 7T = 717/ (2g), population from AO goes 
to Al, and from Bx to Cx, Vx. During the transfer, the 
leakage contributions to the populations in Al and Cx 
thus are first-order in 77-, and so much less than unity. 
Nevertheless, to obtain optimal cooling to first-order in 
"ffi/g we must also optimally transfer these populations 
to the ground state. The population in Al is already be- 
ing optimally rotated to Bl, and since the leakage popu- 
lation is small, and distributed over the great circle path, 
it is clear that the best option is simply to maintain this 
rotation: the presence of thermalization does not change 
the optimal protocol. The populations that appear in 
Bx, Vx, can be rotated to states in the B subspace that 
are not already destinations of existing rotations, if the 
auxiliary is large that additional B-states exist. This tells 
us two things. First, nothing can be done to retrieve the 
populations in B states if the auxiliary is only as large 
as the target. Second, once the Hamiltonian is chosen 
to rotate all the C states to the ground state, all the 
first-order leakage is handled optimally, and no further 
cooling is possible (to first-order). 

If the protocol above is indeed optimal, then given the 
structure of the thermal transition rates in Harmonic os- 
cillators and single qubits, the maximum cooling can be 
achieved with an auxiliary dimension M = 2N — 1, and 
the majority of the cooling with M — N. Further, the 
optimal interaction Hamiltonian, for an auxiliary with 
dimension M, would be H^ pt = G + , where 
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and |m, n) means target state \m) and auxiliary state \n). 

Using the above cooling Hamiltonian we can now cal- 
culate explicitly the final ground-state population for a 
single qubit and a harmonic oscillator. To do this we use 
the linear version of the quantum-jump stochastic mas- 
ter equations that are equivalent to the thermal master 
equations [IB] . This allows us to exploit the fact that all 
the dynamics to first-order in 7T is captured by trajecto- 
ries containing no more than one jump. Further, we find 
that to first-order, the maximum ground-state popula- 
tion is reached at t = r (the fact that the true maximum 
is slightly before r is a second-order effect). For a single 
qubit, the minimum population left in the excited state 
is 
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and Pp is the excited-state population at the ambient 
temperature. Thus if Pt is initially small, the fractional 
reduction in the excited-state population (the approxi- 
mate "cooling factor") is simply 7T7/(4(j). 

For the harmonic oscillator, the minimum population 
that remains outside the ground state is 
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where n is the average number of phonons/photons in the 
target at the ambient temperature. When n is small, the 
ground state population at the ambient temperature is 
1 — n, so the approximate cooling factor is, as one would 
expect, the same as for the two-level system. 

We now turn to the second part of our analysis, that of 
using numerical optimization to find optimal cooling pro- 
tocols (we use the BFGS gradient search method |17j). 
We perform these optimizations with a range of target 
and auxiliary dimensions. This analysis confirms, as far 
as it can, all of our assertions above: i) our protocol is 
optimal; ii) no more than 2N + 1 auxiliary states are 
required for cooling resonators; iii) the best cooling is 
obtained just prior to t = n/(2g). Further, we find that 
up to unimportant transformations, the numerical opti- 
mization always finds our protocol, suggesting that it is 
unique. In Fig.[2]ja) we show results for a two-level target 
with a three-level auxiliary (n — 0.5), and two cases of a 
target and auxiliary both with four levels, with h = 0.5 
and n = 0.1 (The latter gives us results for a resonator.) 
These show that the optimization agrees with our proto- 
col. The inset shows cooling using our protocol up to the 
optimal time (for (N,M) = (2,3)), and the circles give 
numerical optimization past this time, confirming that 
no better cooling is achieved for t > n/(2g). We have 
also run the optimization for larger auxiliaries. All these 
results use 7 = 0.01, g = 1. 



All our analysis so far has assumed an ideal auxiliary 
system, since we are interested in the absolute limit to 
cooling. But nanoresonators are presently cooled via 
auxiliaries with significant damping, so finding optimal 
protocols with damped auxiliaries is an interesting ques- 
tion for future work. In Fig. [2](b) we compare our proto- 
col to numerical optimization, when the auxiliary is very 
strongly damped (damping rate k — g). Remarkably, for 
a two-level system our protocol appears to remain op- 
timal, and the maximal cooling remains essentially the 
same. This is very interesting because this also applies 
to resonators in the low-temperature limit. For larger 
systems, we find that our cooling protocol is no longer 
optimal, and the optimized cooling is also degraded by 
the damping, as we would expect. 
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